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ABSTRACT 

This  thesis  investigates  the  synthesis  of  smoothed  impulse  responses  using  sampled 
data  of  the  input  and  output  of  random  noise  driven  electromagnetic  systems.  Special 
interactive  software  was  developed  for  NPS's  time  domain  electromagnetic  scattering 
laboratory.  The  system  performs  signal  acquisition,  synthesizes  time  and  frequency  do- 
main scattering  responses  using  broad  band  random  noise  and  pro1. 'ides  results  as  easily 
evaluated  graphic  displays.  Attempted  validations  of  the  system  arc  made  by  comparing 
synthesized  impulse  responses  for  microwave  filters  and  transient  scattercrs  to  alternate 
experimental  and  computational  data. 
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I.     INTRODUCTION 

A.      OVERVIEW 

The  objective  of  this  research  is  to  demonstrate  the  viability  of  performing  high- 
resolution  impulse  response  scattering  measurements  using  a  broad-band  noise  source. 
The  development  of  laboratory  facilities  for  high-resolution  impulse  response  scattering 
measurements  has  generally  proceeded  using  two  major  techniques.  The  first  of  these 
employs  a  stepped-frequency  coherent  oscillator  and  a  vector  (magnitude  and  phase) 
receiver.  Impulse  response  target  characteristics  are  obtained  via  inverse  Fourier  trans- 
formation of  the  frequency  domain  data.  The  second  method  obtains  impulse  response 
measurements  directly  in  the  time-domain  by  use  of  a  repetitive  fast-pulse  target  illumi- 
nation with  a  sampling-scope  acting  as  the  receiver.  The  time-domain  approach  for 
scattering  measurements  oilers  a  viable  alternative  to  the  more  prevalent  continuous 
wave  approach.  Transient  scattering  measurements  provide  waveforms  that  can  be  more 
directly  interpreted  as  to  cause  and  effect  and  allow  exact  range-gating  of  target  re- 
sponses for  elimination  of  unwanted  clutter  [lief.  1.  2j. 

The  development  of  digital  sampling  oscilloscopes  and  broad  band  noise  generators 
made  possible  a  third  method  for  high-resolution  impulse  response  scattering  measure- 
ments. This  technique  is  based  on  the  estimation  of  the  crosscorrelation  function  be- 
tween input  and  output  using  a  broad-band  noise  source  as  the  transmitter.  This 
methodology  has  been  employed  in  measuring  the  mechanical  impulse  response  of  large 
structures,  such  as  bridges  and  buildings  [Ref.  3). 

The  practical  advantages  of  this  third  technique  are  two-told.  The  fust  is  the  wide 
bandwidth  and  high  power  available  from  noise  sources  which  are  available  at  much 
lower  cost  than  comparable  stepped-frequency  sources.  A  conclusive  demonstration  of 
the  viability  of  noise  source  impulse  response  measurements  may  lead  to  further  refine- 
ments and  ultimately  to  commercial  marketing  of  this  technology.  The  second  advan- 
tage is  related  to  the  use  of  noise-source  illumination  lor  tactical  and  strategic  radar 
applications.  An  obvious  benefit  would  be  the  masking  of  the  radar  interrogation  signal; 
this  would  appear  at  the  target  as  either  jamming  or  interference.  Confusion  would  re- 
sult as  to  proper  countermeasurcs  to  be  employed  by  the  target.  Furthermore,  since  this 
method  uses  random  noise,  which  is  orthogonal  to  other  signals  and  other  noise 
(producing  zero  crosscorrelation  with  these),  a  high  SXR  requirement  could  be  achieved. 


This  thesis  is  a  continuation  in  a  series  of  efforts  in  transient  electromagnetic  scat- 
tering that  began  in  1979  at  the  Naval  Postgraduate  School  (NPS).  The  development 
of  digital  sampling  techniques  allowed  the  development  of  transient  scattering  ranges  in 
the  late  1960s  [Ref.  4].  A  transient  scattering  range  having  sufficient  bandwidth  and 
signal  to  noise  ratio  (SNR)  to  support  radar  target  identification  was  initially  con- 
structed at  NPS  in  1980  using  a  ground  plane  configuration  [Ref.  5,6].  A  very  wideband 
(short  pulse)  free-field  scattering  range  was  then  constructed  in  1983  for  implementation 
of  target  identification  based  upon  natural  resonances,  as  introduced  by  Mittra  and  Van 
Blaricum  [Ref.  7].  This  free  field  transient  range  was  validated  by  Mariategui,  and 
McDaniel  [Ref.  8,9]. 

B.      NOISE  SOURCE  IMPULSE  RESPONSE  MEASUREMENT  CONCEPT 

In  the  noise  source  impulse  response  measurement  technique,  a  broad  band  noise 
signal  is  used  as  the  input.  The  system  impulse  response  is  derived  from  measuring  the 
crosscorrelation  function  between  the  input  and  output. 

The  configuration  shown  in  Figure  la  lias  been  used  to  implement  the  noise  source 
impulse  response  measurement  technique  in  the  NPS  scattering  laboratory.  It  was 
thought  at  first  that  the  Tektronix  7S04  Digital  Processing  Oscilloscopes  (DPO)  would 
allow  simultaneous  measurements  of  both  channels.  After  much  effort,  this  was  found 
to  be  untrue,  as  will  be  shown.  A  new  method  was  developed,  using  the  measurement 
of  the  difference  of  the  input  and  output  time  sequence.  A  digital  estimator  of  the 
crosscorrelation  function  is  formed  in  the  data  processor  for  the  measurement  and  the 
expectation  of  the  crosscorrelation  function,  R,y[n).  is  obtained  by  further  averaging  oC 
this  estimator  using  N-block  averages  of  the  stored  data. 

The  complete  scattering  measurement  using  a  noise-source  requires  three  steps. 
Looking  at  the  scattering  range  system  representation  in  Figure  lb,  the  target  transfer 
function  II.(J)  will  be  evaluated  by  measuring  X(n)  and  Y(ii),  the  input  and  output  of  the 
respective  transmitting  and  receiving  antennas.  The  second  step  necessitates  factoring 
out  the  effects  of  the  antennas,  in  addition  to  eliminating  the  signal  pollution  of  the 
spurious  cross-coupling  and  clutter  in  the  chamber.  The  following  three  measurements 
are  made  with  equivalent  svstem  diagrams  shown  in  Figure  2: 
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B.  System  Diagram 
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Figure   1.      The  Measurement  of  Noise  Source  Impulse  Response 


Step   1.   :   Background  measurement  (no  target) 
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H2(f)    =    Hi(f)    +   Ht(f)   HP(f)   Hr(f) 
Test  target  measurement 
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Step  4.   :  Subtraction  of  the   background  effect 
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Hs(f)   =   H3(f)   -   Hi(f) 
Step     5.   :   Extraction  of  the  target  response   (Deconvolution) 
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Figure  2.       The  Procedure  of  Noise  Source  Impulse  Response  Scattering  Measure- 
ment 


1.  No  target  present,  measure  /i((r) 

2.  Calibration  target  (sphere)  present,  measure  h2(t) 

3.  Test  target  present,  measure  /z3(f) 

Then,  /?,(/)  will  be  subtracted  from  h2(i)  and  h3{t)  to  eliminate  the  clutter  effect. 

lh{t)  =  }h{t)-h{{t)  (1.1) 

h5{t)  =  h2{t)-h{{t)  (1.2) 

Finally,  the  desired  impulse  response  is  extracted  from  /?4(0  and  /i<U)  using  an  optimal 
deconvolution  technique  known  as  Riad's  method  [Ref.  10]  which  will  be  described  in 
detail  in  the  following  chapter. 

An  important  difficulty  encountered  in  this  thesis  research  was  that  no  sampling 
device  with  a  high  enough  rate  for  the  source  signal,  having  a  bandwidth  of  approxi- 
mately 13  GHz,  is  available  thus  far.  Consequently,  the  estimation  of  the  crosscorre- 
lation  function  could  not  be  computed  using  a  pair  of  properly  sampled  sots  of  the  time 
function  of  each  input  and  output  signal.  Two  alternate  methods  were  developed  which, 
for  the  case  of  ergodic  random  noise,  permits  arbitrarily  slow  sampling  to  be  used.  These 
two  Nyquist  independent  techniques  will  be  described  in  Chapter  II. 

C.      OVERVIEW  OF  THESIS 

The  objectives  of  this  thesis  were  to: 

1.  Develop  a  working  software  program  which  will  acquire  transient  response  data 
from  the  target  and  compute  the  estimation  of  the  crosscorrelation  function  fol- 
lowed by  a  computation  of  the  system  impulse  response  using  a  deconvolution 
technique. 

2.  Demonstrate  impulse  scattering  response  measurements  using  a  broad  band  ran- 
dom noise  source,  and  verify  the  performance  in  comparision  with,  another  tech- 
nique:   time-domain  measurement  using  modified  step  function  input. 

Chapter  II  will  expand  the  theory  of  noise  source  impulse  response  measurement. 
This  will  examine  the  basic  theory  of  analog  crosscorrelation  techniques  lor  acquiring 
system  impulse  responses  and  consider  an  appropriate  derivation  lor  the  discrete  version 
of  the  technique  for  the  sampled  signal.  In  addition,  it  will  quantify  the  estimation  of 
the  crosscorrelation  function. 

Chapter  III  describes  the  experimental  system  in  detail.  '1  his  will  incluo  the  de- 
scription of  the  original  Digital  Processing  Oscilloscope  (DPO)  and  modifications  that 


were  made.  A  description  of  the  noise  source  hardware  and  its  modified  hardware  will 
also  be  given  in  this  chapter.  In  addition,  the  problems  encountered  in  this  research  will 
be  discussed  and  some  guidelines  for  the  laboratory  work  will  be  summarized. 

Chapter  IV  contains  the  calibration  and  validation  measurements  for  simulated 
targets.  This  will  include  the  initial  quality  tests  of  simultaneous  channel  sampling,  noise 
source,  and  crosscorrelation  function  that  were  made  by  sampled  data  using  the  noise 
source.  The  derivation  of  the  crosscorrelation  of  the  measurement  itself 
(autocorrelation)  and  frequency  response  will  be  considered  using  FIT  techniques.  The 
initial  test  involving  a  microwave  bandpass  amplifier  using  a  noise  source  and  a  step 
generator  with  amplifier  will  also  be  discussed. 

Chapter  V  will  describe  the  electromagnetic  scattering  measurement.  The  exper- 
imental setup,  scattering  range  and  characteristics  of  the  antenna  used  will  be  explained. 
Finally,  the  impulse  response  measurement  of  a  metal  sphere  target  will  be  attempted 
using  the  noise  source.  The  resultant  failure  of  this  measurement  will  be  considered. 

Chapter  VI  provides  some  conclusions  about  this  experimental  approach  involving 
noise  source  impulse  response  measurements.  Additionally,  recommendations  are  made 
regarding  improvements  in  the  current  system  and  possibilities  for  further  research. 


II.     THEORY  OF  NOISE  SOURCE  IMPULSE  RESPONSE 

MEASUREMENT 

A.     THEORY  OF  ANALOG  SIGNAL  CROSSCORRELATION 

Looking  at  the  system  diagram  in  Figure  lb.  the  responses  of  transmitter  antenna 
ht(t)  and  receiver  antenna  h,(t)  are  cascaded  with  the  parallel  responses  of  the  target  h.(t) 
and  the  clutter  h£t),  forming  an  overall  response,  h(t), 

h(t)  =  ht(t)  *  [  hs(r)  +  hc{t) }  *  h,{t)  (2.1) 

The  total  measurement  system  can  be  represented  as  a  simple  linear  system  whose  re- 
sponse, \\t),  due  to  an  input,  X(t),  is  given  by  the  convolution 

}'(/)  =  h(t)*X(t) 

(2.2) 
h{o)X{t-  a)  da 


U  can  be  shown  that  the  crosscorrelation  of  the  input  and  output  of  a  linear  system 
estimates  the  system  impulse  response  when  the  input  has  a  bandwidth  that  is  large 
compared  to  the  bandwidth  of  the  system.  This  results  from  superposition  as  applied 
to  stochastic  expectation, 

Rxy(t)  =  ElX(T)Y(T  +  i)\ 

=  E[h(t)*X(r)X(T  +  t)] 

=  h(t)*E[X(T)X(T  +  n) 

Rxy{0  =  h(t)*Rx(t)  (2.3) 

For  ergodic  processes,  the  mean  values  and  moments  can  be  determined  by  lime  aver- 
ages as  well  as  by  ensemble  averages. 

R,xy([)=  <X(z)Y(z  +  n> 

(2.4) 


A'(r))'(r  +  t)dv 
Replacing  the  output  signal  \\t)  by  the  convolution  defined  in  Equation  (2.2)  results  in 
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Figure  3.      Crosscorrelation  Measurement  System  Model 
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(2.5) 


This  leads  to  an  important  property  of  the  Fourier  transform  relation,  which  is 

SxyV)  -  //(/I  Sx(/)  (2.6) 

known  as  the  crosscorrelation  theorem  [Ref.  11]. 

D.     CROSSCORRELATION  MEASUREMENT 

1.     Method  I.  :  Simultaneous  Dual  Channel  Measurement  with  Sampling  Rate 
Greater  than  Nyquist  Frequency. 

Three  methods  of  performing  the  measurement  could  be  used  to  estimate  the 
crosscorrelation  of  the  linear  system.  One  possible,  probably  the  most  general,  way  is 
to  directly  sample  the  input  and  output  sequence  simultaneously  and  store  the  data  in 
the  computer  memory.  This  involves  a  delay  of  the  input  sequence,  and  a  compulation 
of  the  average  of  the  product  of  the  two  sequence  vectors. 


K-l-n 

&xy{n)  =  ~v      l-j       k  ^k+n 
k=n 

K-l-n 
y«)-7    X  X(k)Y(k  +  n)  (2.7) 

k=n 

This  scheme  is  illustated  in  Figure  4  and  compared  to  other  methods.  Two  conditions 
must  be  satisfied  to  use  this  method, 

1.  Sampling  must  be  done  with  a  sampling  rate  greater  than  the  Nyquist  frequency. 

2.  Sampling  of  two  channels   must  be  done  simultaneously. 

2.     Method  2.  :  Simultaneous  Dual  Channel  Measurement  with  Arbitrary  Sampling 
Rate.  (Less  than  Nyquist  Frequency) 

Method  1  could  not  be  used  because  the  required  bandwidth  of  the  DPO 
sampler  is  insufficient.  An  alternative  method  requires  manually  shifting  the  input  (or 
output)  signal,  being  incremented  by  the  sampling  interval.  Samples  of  the  input  and 
output  signal  are  taken  with  a  suitable  sampling  rate,  since  each  sampled  lime  sequence 
is  treated  as  the  sample  set  of  the  ensemble  space  rather  than  a  time  sequence. 

K-] 

*  1    Y"1 

fc=0 
K-\ 

&^(«)  =  ^£aw<)*W<  +  »)  (2.S) 

where, 

n  =  nT  with       y>2/max  (2.") 

TssLT  ,  L$>\  (2.K)) 

flete,  T  denotes  the  time  period  of  delay  and  Ts  is  the  sampling  period.  Using  this 
technique,  the  estimation  of  the  crosscorrelation  of  the  delayed  time  point  can  be  com- 
puted by  each  measurement  of  input  and  output  signal  with  an  appropriate  amount  of 
delay.  This  method  is  illustrated  in  figure  4.  The  remaining  condition  which  must  be 
satislied  for  the  crosscorrelation  measurement  is  simultaneous  sampling. 


For  the  case  of  random  noise,  it  looks  like  the  direct  measurement  method  under 
the  high  rate  sampling  environment.  The  second  method  was  initially  chosen  for  this 
research. 

3.     Method  3.  :  Pre-subtracted  Single  Channel  Measurement  with  Arbitrary 
Sampling  Rate.  (Less  than  Nyquist  Frequency) 

During  the  laboratory  work,  it  was  shown  that  an  unknown  critical  noise  is 
added  to  the  left  sampling  channel  of  the  DPO  which  makes  the  second  method  impos- 
sible to  use.  This  may  be  caused  by  mispositioning  of  the  actual  data  in  the  bulfer.  This 
phenomena  and  the  equivalent  model  are  illustrated  in  Figures  5  and  6.  An  alternative 
method  which  avoids  this  problem  is  to  subtract  Y  from  A'  using  the  built  in  "Add"  and 
"Invert''  functions  in  the  DPO.  This  results  in  only  one  sequence  of  pre-subtracted  data 
and  is  followed  by  computing  the  estimation  of  the  mean-square  value  of  this  data  se- 
quence. 

Referring  to  Figure  6.  assume  that  the  impulse  responses  of  the  sampler  chan- 
nels of  the  DPO  are  unit  impulse  responses.  (The  effect  of  the  sampler  responses  of  the 
DPO  will  be  discussed  in  the  later  section  of  this  chapter.) 

p(t)  =  q{t)=*5(t)  (2.11) 

Then  the  measured  sequences  Xf(n)  and  Yq(n)  can  be  denoted  as 

Xp(n)  s  X(n)  +  Dx{n)  +  Nx(n)  (2.12a) 

Yq{n)s  Y{n)  +  Dy{n)  +  Ny(n)  {2Mb) 

where, 

DJ/i)   :  DC  bias  due  to  the  DPO  vertical  adjustment. 

.Vv(//)     :  other  noise,  (thermal  noise,  jitter,  quantization  noise  etc.) 
with  similar  representations  for  DSn)  and  NSn)  usig  the  F-chaunel. 
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A.  Method  1 

N  samples 

y     MliillMlliilll  .        ,  *-- 

X        I  I  I   I  I  I  I  I  |  I  I  ]  I  I  I   I  RxyW  -  "F     >     A'(/c)  !'(*  +  n) 

A  IT   I    l    !   I   i   I    I    I    I   I    I    II    I  '  A      *— « 

— Hn~^h* —        ~*W*~ 

Delay         Sampling  period 
Dual  channel  measurement  with  sampling  rate  greater  than  Nyquist  frequency. 

B.  Method  2 

N  samples  „  , 

A  — 1 

Y  — \ 1 1 1 R^^-^yXikDYikL  +  n) 

Y  I I I I 

X j — 

T  s  LT 
-*jnTs(—  -*|      T      \+- 

Deiay  Sampling  period 

Dual  channel  measuremient  with  sampling  rate  less  than  Nyquist  frequency. 

C.  Method  3 

N  samples  


Y         I  I  I  |  RjM  =  f  [  [  Z(oo)  f  - 1  Z{n)  f  ]    , 

x    — I 1 1 h- 

1  Z(n)  =  Xp{kL)  -  Yq{kL  +  n) 

— HnTs(—  — )      T      p— 

Delay  Sampling  period 

Single  channel  measurement  with  sampling  rate  less  than  Nyquist  frequency. 


Figure  4.       Comparison  of  The  Three  Estimation  Techniques  for  The  Crosscorre- 
lation 


11 


A.  Subtracted  by  computer 
(after  sampled) 


B.  Subtracted   by  DPO 
(before  sampled) 
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Figure  5.      Simultaneous  Channel  Sampling  Test  by  Subtraction 
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Figure  6.      A  Detailed  Model  of  The  Sampling  System  of  the  DPO 
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Let  a  subtracted  sequence  be  denoted  as 

Z(n)  =  Xp{kL)-  Yq(kL  +  n)  1.2.13) 

where  "n"  is  due  to  the  time-shift.   Then  the  mean-square  of  Z(n)  is  described  as 


[Z(n)\2  =  {Xr(kD-  YjkL  +  u)]1 


=  [{  X(kL)  +  Dx[kL)  +  nx[kL)  }  -  {  Y{kL  +  n)  +  Dy(kL  +  n)  4-  ny{kL  4-  n) 
=  DC  term  -  2Rxy{n) 


=  [Z(oo)f-2Rxy(n) 

where  we  have  assumed  stationarity  of  the  random  process  in  X  and  Y  .    Therefore,  the 
desired  estimation  of  the  crosscorrelation  function  is 


Rxy(n)  =  j  [[Z(oo)Y-  [Z(n)Y\  (114) 

The  results  of  the  trial  measurements  reveal  that  this  estimate  has  a  better  and  more 
consistent  value  than  the  result  from  method  2.  This  measurement  technique  lias  been 
used  for  this  research. 

C.     PROPERTIES  OF  THE  ESTIMATED  CROSSCORRELATION 

Since  the  crosscorrelation  function  has  an  important  role  in  the  analysis  oi  linear 
systems  with  random  inputs,  an  important  practical  problem  is  that  of  determining  the 
quality  of  the  estimation  of  the  function  lor  experimentally  observed  random  processes. 
In  order  to  evaluate  the  quality  of  this  estimate  it  is  necessary  to  determine  the  mean 
and  the  variance  of  Rxy(n),  since  the  estimated  crosscorrelation  is  a  random  variable 
whose  precise  value  depends  upon  the  particular  sample  function  being  used  and  the 
particular  set  of  samples  taken.  The  mean  value  for  Equation  (2.8)  and  (2.14)  can  be 
computed  as  follows: 
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E[Rxy(")]  -  E[  y  YjX^L)  Y{kL  +  n)  ] 
fc-1 
=  -^y£[.Y(A-Lj  Y(kL  +/i)] 

fc=0 


£l4-(/')]  =  ^(/0  (2.15) 

Thus,  the  expected  value  of  the  estimate  is  the  true  value  of  the  crosscorrelation  function 
and  is  an  unbiased  estimate  of  the  crosscorrelation  function.  The  variance  is  denoted 
as 


l'ur[Rxv{n)}  =  E[{Rjn)-Rxv(n)}2} 


E\  {Rxy{n)Y\  -2E\  Rxy{n)  Rxy(n)  \  f  E[  {Rxy(n)Y  \ 

E[[Rxv(n)}2}-{Rxv(n)}2 


=  E[  {  y  S-V(A/-}  Y{kL  +  "]  ]2 !  ~  tRx>M 

k=0 


Var[  Rxy(n)  ]  =  -^  £  ^T  E  [{  X(jL)  YijL  +  n)  }{X{kL)  Y(kL  +  n)  }  ) 

K    J=o  *=o  (2.16) 


"  (  Rxy(n)  f 

Since  j  and  k  are  dilTerent  samples  in  the  ensemble  space,  it  is  reasonable  to  assume  that 
they  are  statistically  independent  random  variables  vvheny'^  k.  Hence.  Equation  (2.16) 
becomes 

E[{  X(jL)Y{jL  +  n)  }  {  X(kL)Y{kL  +/?)}]=  E[{  X{KE)Y(kL  +  n)  }2  |  .  j  =  k 

=  R  i  :(//) 

A    V    v 
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£  [{  X{jL)  Y(jL  +  n)}{  X{kL)  Y(kL  +  n)}}  =  {  E[X(kL)  Y{kL  +  n)\  y  ,j+k 

*  [Rxy(n)}2 
=  B 

where  A  is  the  crosscorrelation  of  X  and  P  while  B  is  the  square  of  the  crosscorrelation 
for  X  and  Y.    Using  this  result,  the  variance  leads  to 

I  ar[  Rxy(n)  j  =  -~  [  KA  +  (  K2  -  K  )  B2  }  -  D2 
A-B2 


1        2 


Var[Rxy(n)]  =  -i?o2xv(n)  (2.17) 


where  a2ry{n)  is  the  true  variance  of  the  product  of  the  random  variables  X  and  Y:.nT. 

An  interesting  phenomena  takes  place  when  the  time  goes  to  infinity.    In  this  case. 
if  X(n)  and  Y(n)  are  zero  mean  random  variables,  then 

/?Jcy(oo)  =  I2F2 
Rxy(oo)  =  0 

Therefore, 

— 2  — 2 

Var{RxAoo)\=^-$-  (2.18) 

A. 

This  means  that  the  variance  of  the  measured  crosscorrelation  approaches  a  constant 
value  as  the  time  separation  goes  to  infinity.  Therefore,  the  resultant  estimated  cross- 
correlation  could  be  thought  of  as 


Rxy{n)  =  Rxy{n)  +  — —  0{n)  1 2. 1 9) 

where.  0(n)  is  the  estimation  error  with  unit  variance.  However,  this  effect  can  be  re- 
duced by  increasing  the  number  of  sample  points,  since  the  variance  is  normalized  by 
this  number. 
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D.     SYSTEM  IDENTIFICATION  USING  ESTIMATED  CROSSCORRELATION 
FUNCTION. 

To  measure  the  input  and  output  noise,  some  kind  of  measuring  device  must  be 
used.  Therefore  it  would  be  reasonable  to  think  that  the  measured  sequence  results  from 
a  convolution  of  the  sequence  to  be  measured  and  the  system  response  of  the  measuring 
device.  This  is  illustrated  in  Figure  7.  Therefore,  the  measured  crosscorrelation  function 
is  Rx ,  (t)  rather  than  RJt). 


Rx ,  (/)  = 


where, 


pOO 


XD{T]}-(v  +  t)dv 


Xp{t)=  I™  p(o)X(x  -  a)  do 

=  p{x)  *A'(r) 


(2.21 


W  = 


q(fx)  Y{x  -  fi)  d\x 


=  q(r)*  Y{x) 


Substituting  Equations  (2.21)  and  (2.22)  into  Equation  (2.20)  yields 


(2.22) 


J»oo 
— o 


p(o)x(x  —  a)  da 


q{/J.]y(x  +t  —  n)  d/.idx 


P(°) 


PKo) 


oo 


oo 


.\(r  —  a)y{x  +i  —  /<)  dx  q(p.)  dfxda 
x(a.)y((t  +  a  —  fi)  +  a)  da  q(/j.)  d/.ida 

—or:  J  — oo 


jO^I  — OO 

oo     />  oo 


Rxy(i)  = 


pin) 


RjV(t  +  a  —  u)tj{/.i)  dfxda 


Since  the  inner  integral  of  the  right  expression  in  Equation  (2.23)  is  the  convolution  be- 
tween R„y{fi)  and  q(fi),  the  equation  can  also  be  denoted  as  follows, 
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Figure  7.      A  Detailed  System  Model  to  be  Measured 


&x  v  (0  =         P(°}         q{t  +  a  -  ix)RXY{n)  d\ida 
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J  —CO 

Therefore,  the  resultant  measured  crosscorrelation  is 

flVf(r)  =  Rpq(t)  *  Rxy(t) 
=  Rpq{t)  *  /?x(/)  *  MO 

=  RxAi)  *  h{t) 


(2.241 


(2.25 


where  Kx ,  (/)  is  referred  to  as  the  system  crosscorrelation  in  this  thesis  and  has  the  same 
role  as  the  autocorrelation  in  Equation  (2.3). 

Since  the  measurement  of  the  system  crosscorrelation  is  performed  without  the  tar- 
get, it  should  be  an  even  function  if  the  responses  of  the  two  sampling  systems  of  the 
DPO  are  exactly  the  same.  In  fact,  the  responses  of  the  two  channels  do  not  have  exactly 
the  same  characteristics  so  the  test  result  of  the  measured  system  crosscorrelation  has 
an  almost  even  shape,  as  expected.  This  will  be  described  in  chapter  4. 

The  Fourier  transform  of  Equation  (2.25)  results  in 

sv,  (/)  =  spq(f)  sxy(j) 

=  Spq(J)Sx(j)lI(J)  (2.26) 

=  Sx  ,  (J)  //(/) 


which  satisfies  the  crosscorrelation  theorem. 

To  get  the  impulse  response,  Ry    (t)  must  be  deconvolved  by  R,  T  (t). 
least  two  measurements  must  be  performed  for  doing  time  decorivolution: 

RXpX(n)  =  Rpq{n)  *  Rx(n) 

Rx„(n)  =  RM*Rx{n)*h(n) 

The  power  spectrums  of  Equations  (2.27,28)  are 


here! ot e  at 

(2.27) 
(2.28) 
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A  A 


SxAk)  =  SpJk)  Sx(k)  (2.29) 


'V"<7 


Sx^(k)  =  Spq(k)Sx(k)H(n)  (2.30) 

Then,  the  frequency  response  of  the  system  is  extracted  by  dividing  Equation  (2.30)  by 
Equation  (2.29). 

5V  v  (A) 
H(k)  =  ~^— 

(2.31) 
SpJk)  Sx(k)  H{k) 


Spq(k)  Sx(k) 

However,  this  frequency  division  must  be  compensated  since  the  estimated  ero^s- 
correhition  function  is  corrupted  by  noise,  as  described  in  Equation  (2.19),  so  that  the 
power  spectral  density  is  biased  by  the  noise  power.  Hie  technique  used  here  for  noise 
compensation  is  known  as  Riad's  Method.  The  "optimal  estimator"  prescribed  by  Iliad's 
Method  [Ref.  10]  is 

Sx.,(k)Sxx{k) 

Hik)  =  —J-L LI (2.32) 

\Sxx{k)\2  +  kWx 

where  Wx  represents  the  average  spectral  power  density  of  the  system  crosscorrelation 
function: 


K-\ 

k-0 


The  parameter  X  is  denoted  as  the  smoothing  parameter  and  establishes  the  lower  limit 
of  the  denominator,  hence  reducing  the  noise  effect  on  the  computed  target  response 
when  the  input  signal  spectrum  is  small.  The  impulse  response  is  derived  by  performing 
an  inverse  Fourier  transform  on  the  target  frequency  response 
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h{n)  =  F-\H{k)}  (2.34) 

Another  important  point  to  keep  in  mind  is  that  we  use  an  FFT  algorithm  which 
has  quite  different  properties  from  the  continuous  Fourier  transform.  Ordinarily,  when 
the  FFT  is  used  to  estimate  the  frequency  response,  windowing  is  performed  to  smooth 
the  result.  Howrever,  applying  a  window  interferes  with  the  circular  convolution  prop- 
erty. Therefore,  in  this  research,  we  will  not  apply  windowing  before  performing  the 
FFT  algorithm.  Instead  of  windowing,  the  crosscorrelation  function  R, ,  (n)  will  be  "zero 
padded"  to  avoid  the  wrap-around  problem  which  comes  from  the  circular  cons'olution 
property. 
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III.     BASIC  EXPERIMENTAL  CONFIGURATION 

A.     DESCRIPTION  OF  DPO  AND  ITS  MODIFICATION 

To  meet  the  goals  of  this  thesis  work,  an  ultra  wide  band  noise  source  and  a  dual 
channel  sampling  device  with  high  sampling  rate  and  precise  synchronization  must  be 
used.  A  dual  channel  synchronous  sampling  Digital  Processing  Oscilloscope  (DPO)  and 
a  broadband  noise  signal  generator  were  used  for  this  purpose. 

The  DPO  is  not  an  ideal  sampling  device  since  it  does  not  sample  at  a  rate  exceeding 
the  Nyquist  frequency.  Even  though  the  sampling  rate  of  the  DPO  is  measured  in  MHz, 
it  can  display  a  periodic  repetitive  signal  having  up  to  12  GHz  in  bandwidth  by  rear- 
ranging the  actual  low  rate  sampled  data  using  a  display  buffer.  Hence  the  signal  dis- 
played on  the  oscilloscope  looks  like  it  has  been  sampled  at  over  the  Nyquist  rate.  The 
sampling  technique  of  the  Tektronix  DPO  sampler  S-6  is  illustrated  in  Figure  8.  The 
Tektronix  S-52  pulse  generator  produces  a  fast  repetitive  pulse  train.  A  time  limited 
output  of  a  lumped  network  can  be  displayed  at.  an  apparent  sampling  nue  exceeding  the 
Nyquist  rate  by  using  the  S-52  to  drive  the  network.  This  is  shown  in  Figure  9. 

However,  this  property  could  not  be  extended  to  a  non-repetitive  signal  since  the 
technique  does  not  reconstruct  the  original  shape  of  the  actual  sequence  by  rearranging 
the  sampled  data  after  it  is  taken  by  the  sampling  head.  Consequently,  we  could  not 
utilize  the  virtual  characteristics  of  repetitive  sampling  of  the  DPO  for  this  research.  The 
broad  b?nd  noise  generator  has  about  13  GHz  bandwidth  which  is  far  beyond  the  single 
shot  sampling  rate  of  the  DPO. 

Using  the  measurement  method  2.  as  described  in  Section  II. B,  one  of  the  signals  to 
be  measured  is  manually  shifted  through  the  required  sampling  interval.  The  measure- 
ment process  treats  the  data  as  samples  of  an  ensemble  space  rather  than  a  time  se- 
quence since  our  intention  is  to  compute  the  value  of  the  crossconelation  of  the  shifted 
lime  difference.  Under  this  condition,  there  is  no  restriction  of  the  sampling  rate  and  the 
DPO  can  be  used.  There  is,  however,  the  need  ol  an  additional  variable  delay  function 
on  one  channel  to  make  the  required  time  shift.  This  function  is  carried  out  using  the 
"Delay''  knob  on  the  channel  0  sampler. 
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Figure  8.      The  Sampling  Scheme  of  the  Tectronlx  7854  DPO 
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A.  A  portion  of  a  rectangular  pulse  produced  by  the  S-52  pulse  generator 
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B.  Modified  pulse  filtered  by  the  DMR005  bandpass  amplifier 
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Figure  9.      A  Tulse  Signal  and  Its  Convolved  Output  by  a  Linear  Sy.s..  m 
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This  delay  knob  has  a  10  nanosecond  delay  range  which  is  enough  for  measuring  the 
crosscorrelation  function  for  the  small  targets  that  are  used.  Assuming  that  the  band- 
width of  the  sampled  signal  is  less  than  12  GHz,  the  required  sampling  interval  is 

r5  =  4-  =  40   psec  (3.1) 

Js 

where,  according  to  the  sampling  theorem, 

fs  =  2B=  25   GHz  (3.2) 

To  obtain  this  type  of  resolution  on  the  delay,  a  10  turn  potentiometer  was  used  instead 
of  the  built-in  one  turn  potentiometer  on  the  10  nanosecond  delay  knob.  Tins 
potentiometer  has  a  five  thousand  division  scale  on  its  lace  within  the  full  range  of  its 
delay.  This  modification  is  illustrated  in  Figure  10. 

During  the  laboratory  measurements,  it  was  found  that  the  actual  time  delay  was 
not  exactly  linearly  proportional  to  the  delay  reading.  A  potentiometer  reading  table 
was  then  used  to  mark  the  value  of  the  reading  corresponding  to  each  rime  delay.  A 
sample  potentiometer  reading  table  is  shown  in  Figure  11.  Although  the  actual  delay  is 
dependent  upon  the  potentiometer  reading  table  used  with  each  measurement,  a  reading 
difference  of  about  3.0  represents  approximately  a  40  picosecond  delay. 

To  create  a  potentiometer  reading  table,  an  additional  measurement  using  the  pulse 
source  was  performed  prior  to  the  noise-source  measurement.  After  creating  several  such 
tables,  it  was  revealed  that  the  potentiometer  reading  has  about  +  3  picosecond  reading 
error  due  to  misreading  and  the  line-width  of  the  displayed  pulse  when  it  is  performed 
on  20  picosecond  per  division  time  scale.  This  is  illustrated  in  Figure  12.  One  factor 
which  increases  this  line-width  is  the  jitter  noise.  This  reading  error  increases  as  the  time 
scale  increases  and  makes  the  potentiometer  reading  table  next  to  useless.  Therefore,  a 
linear  increment  of  the  potentiometer  by  1..5  is  usually  acceptable  if  the  total  time  win- 
dow to  be  measured  is  over  200  picoseconds.  In  such  a  case  the  measurement  can  be 
performed  using  a  20  picosecond  per  division  time  scale. 

Another  problem  encountered  during  the  laboratory  work  is  a  scaling  difference  of 
the  sampling  interval  of  the  two  channels  as  well  as  the  drift  of  the  sampling  frame.  A 
comparision  of  (.V—  }")  on  the  20  psec./div.  and  10  nsec./div.  is  shown  in  Figure  13. 

Looking  at  Figure  13,  we  see  that  the  upper  part  has  noise  with  a  constant  envelope 
while  the  other  demonstrates  a  noise  with  a  consistently  growing  envelope.    Therefore, 
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it  would  be  reasonable  to  think  that  the  source  of  the  noise  in  the  upper  figure  conies 
from  the  jitter  noise  of  the  two  channels  due  to  the  random  sampling  time  error  while 
the  other  includes  the  jitter  noise  and  a  scaling  difference  of  the  sampling  interval.  One 
possible  source  of  the  scaling  problem  might  be  come  from  machine  accuracy.  So,  if  the 
two  signals  are  aligned  at  one  position  using  the  potentiometer  delay  knob,  so  that  the 
subtracted  signal  at  the  point  shows  a  minimum  envelope,  then  the  error  caused  by  dif- 
ferent sampling  rates  grows  with  increasing  the  time  differences  from  the  aligned  point, 
producing  larger  error  in  the  envelope.  Therefore,  it  would  be  better  to  sample  the  data 
on  a  20  picosecond  per  division  (nominal  value)  to  minimize  problems  introduced  by  the 
difference  of  sampling  rates. 

During  the  laboratory  wor-k  it  was  found  that  the  images  of  a  pulse  signal  as  dis- 
played on  each  of  the  two  channels  was  gradually  drifting  from  left  to  right  with  differing 
velocity.  This  drift  problem  might  come  from  the  thermal  inconsistency  of  the  DPO 
since  it  can  be  reduced  by  warming  up  for  about  one  hour.  However,  the  elfect  still  re- 
mains. This  problem  requires  completion  of  the  experiment  in  as  short  a  time  as  possi- 
ble. 

B.      NOISE  SOURCE  GENERATOR  HARDWARE  AND  ITS  DEVELOPMENT 

A  model  DMR-005  broadband  noise  generator  was  used  for  this  research.  This 
noise  source  includes  a  noise  generator  and  two  cascaded  broadband  amplifiers,  each  of 
which  has  about  20  dB  gain  with  a  13  GHz  bandwidth.  The  system  output  is  approxi- 
mately 400  mV  in  rms  value  and  the  voltage  histogram  resembles  a  Gaussian  distrib- 
ution. The  block  diagram,  the  frequency  characteristic,  and  a  sample  recorded  output 
of  the  noise  source  are  illustrated  in  Figures  14.  15,  and  16. 

In  the  laboratory  it  was  found  that  a  frequent  change  of  the  physical  environment 
may  lead  to  electrical  inconsistency  resulting  in  output  fluctuations.  A  little  after  he- 
ginning  this  research,  one  of  the  two  amplifiers  used  in  the  noise  source  failed  and  it  was 
sent  to  the  manufacturer  to  fix  the  malfunction.  After  the  amplifier  was  returned  to 
NTS,  an  aluminum  support  was  attached  to  the  noise  source  to  avoid  further  problems. 
The  support  also  acts  as  a  heat  sink.  Thereafter,  a  consistent  noise  output  has  been 
produced. 
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Figure  10.      Modified  Tektronix  7854  DPO 
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Figure   11.      Potentiometer  Reading  Table 
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Figure   12.      The  Reading  Error  Due  to  the  Line-width  of  Ramp  Signal 
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A.  (X  —  10  displayed  on  the  20  psec./div.  window 


B.  (X  —  Y)  displayed  on  the  10  nsec./div.  window 


Figure   13.      The  Effect  of  the  Difference  of  Sampling  Interval 
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Figure   14.      Noise  Generator 
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Figure   15.      Frequency  Response  of  an  Amplifier  of  the  Noise  Generator 
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Figure  16. 


A  Recorded  Noise  Signal 
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C.      SUMMARY  OF  THE  PROBLEMS  AND  SOLUTIONS  FOR  THE 
CROSSCORRELATION  MEASUREMENT 

The  problems  for  the  crosscorrelation  measurement  and  the  solutions  described  so 
far  are  summarized  as  follows. 

1.  To  produce  a  noise  source  with  consistent  characteristics,  an  aluminum  heat  sink 
support  was  attached  to  the  noise  source. 

2.  A  sampling  rate  lower  than  the  Xyquist  frequency  and  mispositioning  of  the  noise 
added  to  the  channel  0  (output  measurement  channel)  forced  use  of  the  measure- 
ment method  3. 

3.  The  scaling  difference  of  the  sampling  period  of  the  two  channels  leads  to  sampling 
the  data  on  the  20  picosecond  per  division  time  scale  to  minimize  the  effect. 

4.  The  drift  of  sampling  frame  requires  a  warm  up  period  of  at  least  one  hour  to  reach 
thermal  steady  state  and  then  to  perform  the  measurement  as  last  as  possible. 

5.  To  perform  the  fastest  measurement,  a  1.5  reading  difference  of  the  potentiometer 
is  linearly  increased  regardless  of  potentiometer  reading  table  since  the  experiment 
reveals  that  a  1.5  leading  difference  of  the  potentiometer  makes  a  20  picosecond 
delay  with  affordable  error  which  can  be  compromised  with  machine  precision  and 
visual  reading  error. 

6.  Also,  it  is  recommended  to  perform  the  experiment  as  quickly  as  possible  to  reduce 
the  data  transfer  time,  which  is  proportional  to  the  number  ol  the  samples.  (It  takes 
about  8.5  second  tor  transferring  a  set  of  1024  samples  of  data.) 

Following  these  guidelines,  about  half  an  hour  of  measurement  time  is  needed  for  64 
point  by  2  sets  oC  1024  samples  with  an  additional  one  hour  for  warming  up  the  system, 
which  will  be  the  best  condition  for  the  crosscorrelation  measurement  under  the  current 
laboratory  environment.  However,  the  drift  problem  still  remains  since  at  least  one-half 
hour  is  needed  to  perform  the  experiment. 

Finally,  detailed  tests  revealed  that  the  system  crosscorrelation  at  the  peak  point  has 
the  value  of  about  twenty  times  the  peak-to-peak  value  of  the  error  noise.  This  is  illus- 
trated in  Figure  17.  Therefore,  if  the  value  of  the  original  crosscorrelation  drops  down 
significantly,  then  the  error  noise  will  dominate  the  measured  signal  and  will  make  the 
crosscorrelation  method  impossible  to  use. 
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Figure   17.      TheAmountoftheMeasuredErroriNoi.se 
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IV.     CALIBRATION  AND  VALIDATION  MEASUREMENT 

A.  INITIAL  TEST  OF  SIMULTANEOUS  CHANNEL  SAMPLING 

The  most  important  and  critical  role  of  the  DPO  in  this  research  is  in  simultaneously 
sampling  a  pair  of  time  functions  which  are  probed  at  the  input  and  the  output  port  of 
a  linear  system.  A  useful  validation  test  for  the  simultaneous  sampling  characteristic  is 
to  subtract  the  signal  from  one  channel  (X :  input  sequence  in  this  thesis)  from  the  other 
channel  ( Y  :  output  sequence).  This  will  be  displayed  on  the  DPO  screen  by  using  the 
built  in  "Add''  and  "Inverse"  functions  of  the  DPO.  Using  the  manual  delay  knob,  the 
image  of  the  subtracted  signal  can  be  minimized  using  the  simultaneous  sampling  posi- 
tion. The  potentiometer  reading  at  this  point  (it  must  be  the  n  =  0  point  of  the  system 
crosscorrelation,  Rr  .  (/?))  has  a  value  of  about  200  indicating  that  the  sampling  channel 
0  (output  channel)  has  a  time  delay  of  about  2.7  nanoseconds  less  than  that  of  the 
channel  1.  This  is  illustrated  in  Figure  IS  with  a  comparison  with  the  original  input  se- 
quence. 

Looking  at  the  Figure  IS,  the  left  part  shows  a  noise  with  a  consistent  envelope  ol~ 
small  variance.  The  source  of  the  noise  on  the  left  figure  comes  from  the  jitter  noise  of 
the  two  channels  produced  by  random  sampling  time  errors.  Assuming  that  the  envel- 
ope of  the  jitter  noise  is  proportional  to  the  time  derivative  of  the  source  signal  itself,  the 
noise  source  signal  must  have  more  variance  of  its  jitter  noise  vis-a-vis  the  ramp  source 
signal.  This  results  since  the  time  derivative  of  the  ramp  signal  produced  by  DPO  is 
much  smaller  than  that  of  the  noise  source  indicating  that  the  noise  signal  is  more  vul- 
nerable to  the  jitter  noise  than  the  ramp  signal.    This  property  is  verified  in  Figure  19. 

B.  CROSSCORRELATION  OF  THE  DPO  SYSTEM. 

As  described  in  Section  Il.D.  the  crosscorrelation  between  input  and  output  is  con- 
volved by  the  crosscorrelation  between  the  impulse  response  of  the  two  sampling  heads. 
The  discrete  version  of  eq.(2.25)  is 

=  Rpq{n)  *  Rx{n)  *  h{n)  M.I) 

=  Rxx  (/i)  *  h{/i) 
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A.  Channel  1  -  Channel  0 


B.  Channel  1 
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Figure  18.      Simultaneous  Channel  Sampling  Characteristics  by  Subtraction 
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A.  Jitter  noise  of  the  noise  signal 


B.  Jitter  noise  of  the  pulse  signal 


Figure   19.      Jitter  Noise  of  the  Noise  and  the  Ramp  Signal  by  Subtraction 
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where, 

RVq(n)  m  Rpq{n)  *  Rx(n)  (4.2) 

Here  Rr ,  (n)  is  the  estimation  of  the  svstem  crosscorrelation  and  is  convolved  with  the 
crosscorrelation  of  the  sampling  head  and  the  autocorrelation  of  source  noise  itself.  To 
approximately  measure  the  system  crosscorrelation,  we  simply  feed  the  same  input  signal 
to  each  sampling  head,  as  illustrated  in  Figure  20. 

In  the  ideal  case,  each  sampler  system  impulse  response.  p(i)  and  q{t),  must  have  the 
same  shape,  and  at  best  the  unit  impulse  response.  Although  these  two  impulse  re- 
sponses are  not  exactly  identical,  they  will  in  reality  be  very  close  to  each  other.  Conse- 
quently,  the  shape  of  the  svstem  crosscorrelaiion  R, ,  (n)  is  close  to,  but  not  exactlv.  an 

i  r  .  xpxq\ 

A 

even  function,  while  the  autocorrelation  of  the  source  noise  Rx(n)  must  be  even. 

The  initial  measurements  were  performed  using  a  2.5  picosecond  sampling  period  to 
get  a  more  detailed  time  shape  of  the  system  crosscorrelation.  Subsequent  measurements 
were  performed  using  a  10  picosecond  sampling  interval  to  see  the  entire  shape  of  the 
system  crosscorrelation.  Figure  21  and  22  show  the  test  results.  As  expected,  the  shape 
of  the  crosscorrelation  is  almost  even,  which  means  the  two  sampling  channels  have 
similar  characteristics. 

C.     DECONVOLUTIGN  VALIDATION  TEST 

1.     Test  with  Noise  Source 

An  Avantek  model  SAS3-2954  solid-state  microwave  amplifier  was  used  during 
this  validation  test.  This  amplifier  has  an  average  gain  of  42.5  dB  over  the  band  from  2 
to  6  GHz  and  provides  a  3  dB  bandwidth  of  about  5  GHz.  The  spectral  characteristic 
of  the  Avantek  amplifier  is  shown  in  Figure  23.  "I  his  amplifier  is  linked  to  two  20  dB 
attenuators,  since  it  produces  moderate  power,  to  prevent  any  possible  damage  to  the 
sampling  heads  which  have  an  operating  range  of  5  volts  peak.  Figure  24  illustrates  the 
basic  setup  for  this  test. 

To  extract  the  impulse  response  of  the  Avantek  amplifier  (2-port),  the  following 
three  steps  were  followed: 

1.  Test  without  the  2-port,  measure  the  system  crosscorrelation. 

2.  Test  with  the  2-port,  measure  the  crosscorrelation  of  the  2-port. 

3.  Perform  deconvolution  computation. 
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The  system  crosscorrelation  measurement  was  performed  using  6-4  time  points 
having  a  40  picosecond  increment.  This  gives  about  a  2.5  nanosecond  window,  which  is 
almost  three  times  longer  than  the  effective  period  of  the  system  crosscorrelation.  The 
amount  of  the  delay  time  interval  (40  picoseconds)  is  sufficiently  small  in  this  case  since 
this  gives  a  Xyquist  frequency  exceeding  the  bandpass  frequencies  of  the  Avantek  am- 
plifier. A  plot  of  the  measured  system  crosscorrelation  is  shown  in  Figure  25. 

The  bandpass  amplifier  crosscorrelation  measurement  was  performed  using  the 
same  time  points  and  delay  interval.  Further,  an  initial  delay  was  given  to  the  channel 
0  in  the  amount  of  400  picoseconds  to  compensate  for  the  effective  delay  through  the 
Avantek  amplifier.   The  resultant  measured  crosscorrelation  is  shown  in  Figure  26. 

Finally,  the  computation  of  deconvolution  was  performed  using  Iliad's  optimal 
compensation  technique.  As  described  in  Section  II.D,  this  deconvolution  method 
serves  to  limit  the  lower  bound  of  the  denominator  and  thus  reduces  the  noise  effect. 
Several  computations  were  performed  with  different  smoothing  factors  and  finally  the 
case  of  0.1  was  selected  for  the  best  solution.  After  the  computation  of  deconvolution, 
the  computed  impulse  response  was  moved  32  points  to  the  right  to  show  the  best 
graphic  interpretation.  The  result  is  illustrated  in  Figure  27. 
2.     Test  with  Modified  Pulse  Source 

The  same  three  measurement  steps  were  performed  using  the  time  domain 
method  described  by  McDaniel  [Ref.  9j.  The  output  ramp  signal  from  the  S-52  pulse 
generator  was  filtered  by  one  of  the  amplifiers  used  with  the  noise  generator  to  produce 
a  roughly  similar  power  spectral  density  at  the  input.  This  source  signal  has  enough 
innate  zero  pacding  to  eliminate  the  wrap-around  problem.  The  test  was  fairly  straight 
forward  and  the  results  are  shown  in  Figures  28,  29,  and  30. 

The  comparison  of  the  two  methods  is  given  in  Figure  31.  These  two  results 
have  almost  the  same  shape,  as  is  expected,  aside  from  the  time  scaling  difference.  "1  he 
crosscorrelation  method  uses  about  a  40  picosecond  time  interval  while  the  time  domain 
method  samples  with  a  period  of  39  picosecond,  which  is  fixed  by  the  DPO  spaces. 
Another  fact  is  that  the  result  from  the  crosscorrelation  method  has  more  error  noise 
power  than  that  of  the  pulse  method.  Therefore,  the  performance  of  the  crosscorrelation 
method  using  the  current  DPO  appears  to  produce  lower  quality  results  than  that  of  the 
pulse  method. 
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Figure  20.      DPO  System  Crosscorrelation  Measurement  Setup 


41 


W\ 


r— " 


TTT 


r 


Points  :    128    points 

Resolution    :    2.5   picosecond 
Window  :    320   picosecond 


Figure  21.      DPO  Svstem  Crosscorrelation  Function 
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DPO  System  Crosscorrelation  Function  and  Power  Spectral  Density 
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figure  23.       Frequency  Response  of  the  Avantek  Amplifier 


44 


xft) 


:h1  sampler 

D(t) 


xp(KU 


Noise  generator 


inear  system 
h(t) 


Potentiometer  delay 


<Vt 


y(t) 


(n) 


-*«      n       - 


;h0  sampler 
q(t) 


Trigger  (kU 


y'kL+n) 


+     K- 


xn(kL)-yn(kL+n) 


Computer 


I   R    (n) 


Figure  24.      Noise  Deconvolution  Validation  System  Setup 
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Figure  25.      The  System  Crosscorrelation,  Crosscorrelation  Method 
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Figure  26.      The  Crosscorrelation  of  Avantek  Amplifier,  Crosscorrelation  Method 
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Figure  27.      The  System  Response  of  Avantek  Amplifier,  Crosscorrelation  Method 
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Figure  28. 


The  Modified  Pulse  Input  Signal,  Time  Domain  Method 


49 


KUHBER  OF  SAMPLE  POINT  :  256 


WAUEFQRM 

SUP 

:  256  PTS. 

RES  ' 

:  3?. 059  PS. 

win  : 

;  9999. 104  PS. 

mx   : 

I  8.643  ULT 

hih-  : 

-U.775  ULT 

DYN  : 

1.417  ULT 

BIAS: 

3.868  ULT 

AUG  : 

188  TIMES 

dZtO 


PECi 

"RUM 

SMP 

:  129  PTS. 

RES 

:  9. 1  GHZ 

WIN 

:  12.9  GriZ 

MAX 

:  6. OB?   ULT 

F;gure  29.      The  Output  Signal  of  the  Avantek.  Amplifier,  Time  Domain  Method 


50 


WJMBER  OF  SAMPLE  POINT  :  256 

WAVEFORM 


TEST  WITH 
TARGET 


TEST  WITHOUT 
ThRGET 


TARGET 
RESPONSE 


SPECTRUM 


SMP. 

I  256  PTS. 

SMP. 

:  129  PTS. 

RES.  ! 

:  39.859  PSEC. 

RES. 

:  0.1  GHZ. 

win.  : 

!  9595. 104  PSEC. 

win.  : 

I  12.9  GHZ. 

MAX.  : 

:  6.318  ULT. 

MAX.  . 

;  2.SSS  ULT. 

MIN.  : 

:  -8.387  ULT. 

DVN.  ' 

Q.704  ULT. 

SMGGTf 

4  FCT.  :  u. 1 

Figure  30.      The  System  Response  of  Avantek  Amplifier,  Time  Domain  Method 


51 


IMPULSE  RESPONSE 


♦    :    NOISE 
o   ;    PULSE 


/\nj  AA 


y  ■  r»  4t      lag 


v    v  v 


SAMPLES     :  64  POINT. 
RESOLUTION  :  40  PSEC, 
WINDOW      :  2560  PSEC. 


Figure  31.      The  Comparison  of  the  Tno  Methods 
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V.     ELECTROMAGNETIC  SCATTERING  MEASUREMENT 

A.     DESCRIPTION  OF  THE  SCATTERING  RANGE 

To  derive  the  characteristics  of  the  noise  source  impulse  scattering  response  meas- 
urement, a  physical  environment  and  equivalent  model  for  the  scattering  range  should 
be  developed.  Since  the  scattering  range  was  already  described  in  detail  by  Mariatequi 
[Ref.    8  ],  and  McDaniel  [Ref.  9],  only  a  brief  description  will  be  given  in  this  section. 

The  general  physical  environment  of  the  scattering  range  is  illustrated  in  Figure  la. 
As  shown  in  the  illustration,  the  system  can  be  divided  into  5  parts. 

1.  Anechoic  chamber 

2.  Transmitting  and  receiving  antennas 

3.  Broadband  analog  noise  generator 

-4.   Dual  channel  sampling  device  (DPO) 

5.   Signal  processor  and  controller  (Tektronix  4052A  desktop  compute] ) 

The  purpose  of  the  anechoic  chamber  is  to  shield  the  incident  and  scattered 
waveforms  from  external  noise  and  interference  effects  as  well  as  multipath  contam- 
ination.  The  physical  dimensions  of  the  chamber  are 

1.  Longitudinal  length  :  20  ft. 

2.  Lateral  width  :  10  ft. 

3.  Height  :  10  ft. 

To  provide  shielding  from  atmospheric  and  man-made  noise,  the  chamber  is 
internally  covered  with  aluminum  panels  which  are  earth  grounded.  A  special  absorber 
material  manufactured  by  Rantec,  a  division  of  Lmerson  Electric,  is  attached  to  the 
aluminum  shield  to  absorb  the  electromagnetic  radiation.  This  absorber  is  ridged  along 
the  length  of  the  chamber  to  guide  the  incident  wave  to  the  back  wall  where  the  most 
of  the  incident  and  guided  electromagnetic  wave  is  absorbed  by  18  inch  pyramidal  cones. 
Eight-inch  pyramidal  cones  are  attached  to  the  front  wall,  where  the  antennas  are 
mounted.  These  absorb  back-scattered  radiation  and  prevent  reradiation  of  the 
electromagnetic  wave.  This  material  is  designed  to  provide  back-scattering  attenuation 
of  a  500  MHz  signal  by  about  3«.)  dli  below  that  of  a  (lat  metal  plane.  Reflection  coeffi- 
cients are  increased  to  about  -12  dB  at  1  GHz  for  near  grazing  incidence  on  the  side- 
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walls,  floor  and  ceiling.  This  leads  to  a  realistic  low  frequency  limitation  of  about  1  GHz 
due  to  the  resultant  multi-path  interference. 

A  styrofoam  column  is  used  as  the  support  for  the  various  targets  in  the  chamber. 
This  material  has  a  relative  permittivity  of  about  1.1  at  3  GHz  so  that  it  gives  negligible 
effect  (reflection,  refraction)  to  the  transmitted  wave  but  provides  a  stable  support  for 
the  targets. 

Two  double-ridged  horn  type  antennas  are  used  for  the  transmitting  and  receiving 
antennas.  They  have  a  usable  bandwidth  of  1  to  12.4  GHz  and  have  a  relatively  Hat  gain 
over  the  band. 

The  broadband  analog  noise  generator,  DPO,  and  the  Tektronix  4052A  desk,  top 
computer  were  already  described  in  Chapter  III. 

B.     DERIVATION  OF  NOISE  SOURCE  IMPULSE  SCATrERRING  RESPONSE 
MEASUREMENT  FORMULA 

Since  the  goal  of  (.his  research  is  to  validate  the  crosscorrelation  measurement  by 
comparing  [he  impulse  scattering  response  result  with  that  from  the  time  domain  meas- 
urement, a  simple  sphere  was  chosen  for  the  target. 

An  equivalent  system  diagram  for  the  scattering  range  is  shown  in  Figure  32.  Here. 
the  impulse  responses  are 

1.  h5{i)       :  Target 

2.  he(t)       :  Clutter  and  antenna  cross-coupling 

3.  /;.(/•)       :  Transmit  antenna 

4.  hr(i)       ■  Receiving  antenna 

and  the  signals   are 

1.  x(t)        :  Transmitted  (source)  signal  produced  by  DPO 

2.  y(l)        :  Received  (back  scattered)  signal 

3.  x„(kL)   :  Sampled  transmitted  sequence 

4.  yq{kL  f //):  Sampled  received  sequence 

5.  ziti)         :  Stored  sequence  (  A'-  }') 

6.  Rx    (n)   :  Crosscorrelation  output 

Compared  with  the  discussion  in  Section  IV. C.  the  system  impulse  response  in  the 
anechoic  chamber  is  more  complicated  than  a  simple  2-port  network  response.  The 
impulse  response  of  the  anechoic  chamber  is 
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Figure  32.      A  System  Diagram  of  the  Scattering  Range 
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h(t)  =  hit)  *  (hs(t)  +  hc{i))  *  hr{t) 

=  (A/(/)*Ar(/))*(A,(r)  +  /ic(r))  (5.1) 

=  ha(t)  *  (/MO  +  hc{i)) 

This  equation  suggests  the  need  of  at  least  three  sets  of  measurements  for  extracting  the 
target  response  hs{t).  Using  our  method,  the  three  crosscorrelations  that  must  be  meas- 
ured are 

Rx^{n)x  =  Rpq{t)  *  Rxy{t) 

=  [  Rpq(t)  *  Rx(i)  *  hz{i)  *  hr{t)  ]  *  hc([)  (5.2) 

=  'MO  *  fcc(0 

V«)2  =  ^)MV')  +  /',(')1  (5.3) 

^('Oa  =  0(0  *  I  ^(0  +  W)  1  (5.4) 

where. 

0(0  =  Rpq(t)*Rx(t)*  hfc)*hk)  (5.5) 

However,  experiments  revealed  the  need  of  one  more  background  measurement  for  bet- 
ter subtraction  of  the  background  noise.  This  is  because  of  the  use  of  a  relatively  shoit 
time  window  compared  with  the  difference  of  the  propagation  delays  of  different  targets, 
(the  calibration  target,  and  real  target).  Therefore,  instead  of  Equation  (5.2)  to  (5.4), 
we  use  Equation  (5.6)  to  (5.9). 

Rx    (n)n  =  4>{n)  *  hc{n)2  (5.6) 

^.v./">i3  =  </>(")  *hjin)2  (5.7) 

Rxy  (n)2  =  4>{n)  *  [  hp{n)  +  /!.('?);  J  1 5.8) 


Rxv{n)2  =  <fi(n)  *  [  hs{n)  +  hc(n)3  \  (5.9] 

where 
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<$>{n)  m  Rpq{n)  •  Rx(n)  *  h£n)  *  /ir(/i)  (5.10) 

hc(n)3^hc(n-m)2  (5.11) 

Here,  m  is  the  time  delay  due  to  the  different  measurement  time  origin  caused  by  tiie 
physical  shape  of  each  target.  So,  the  best  delay  could  be  determined  by  the  difference 
of  the  transit  time  of  the  leading  edge  of  the  scattered  electromagnetic  wave  lor  each 
target. 

in  =  Intecer  {  — rr-  } 
cls 

In  this  equation,  D  denotes  the  radial  difference  of  the  two  spheres  and  c  represents  the 
free  space  velocity  of  an  electromagnetic  wave.  This  is  illustrated  in  Figure  33  and  an 
example  of  this  effect  is  shown  in  Figure  3-1  using  an  I  sq.  ft.  copper  sheet  target  for 
better  reflection. 

The  next  tiling  to  do  after  the  lour  measurements  of  the  crosscorrelation  is  the  ex- 
traction of  the  target  response  hs{n).   This  is  performed  in  two  steps. 

1.  Clutter  subtraction. 

2.  Crosscorrelation  deconvolution  using  Riad's  method. 

At  first,  subtraction  of  the  clutter  effect  was  performed.    Let  the  subtracted  result  be 
denoted  as 

A  A  A 

(5.12) 

A  A 

=  4>{n)  *  hp{n) 

R3X(n)  =  R     (n)3-R     (n)u 

(5.13) 

The  following  is  the  deconvolution  for  extracting  the  response  of  the  target.  The 
key  point  of  the  deconvolution  is  the  use  of  the  calibration  target  which  has  a  known 
computed  transfer  function.  The  Mie-serics  program  prepared  by  Professor  Morgan 
computes  the  transfer  function  of  the  sphere  [Ref.  9J. 
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Figure  33.      Path  Differences  for  Target  and  Calibration  Spheres 
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Figure  34.      Examples  of  the  Delayed  Signal  Due  to  Different  Paths 
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Let  HpciJ)  represent  the  computed  transfer  function  of  the  8  inch  sphere.    For  an  ideal 
case,  the  transfer  function  of  the  target  could  be  computed  as 


S31W 


HM  =  -frx  hpM  (5AA) 


where  S31(/)  and  S2l(J)  are  the  Fourier  transforms  of  the  crosscorrelations  R3l(t)  and 
R2l(t).  Since  we  are  dealing  with  the  estimated  crosscorrelations.  Riad's  method,  as  de- 
scribed in  Section  II.D,  is  applied. 

Hs(k)  =  — -± U IIp,(k)  (5.15) 

\S2l(k)\2  +  AW 


21 


where. 


K- ' 

1  Y  ■ 

,r-i=Y  /  s'-'lk)  |V1'" 

fc=0 

Finally,  the  smoothed  estimation  of  the  impulse  response  will  be  obtained  from  the  in- 
verse Fourier  transform  of  Equation  (5.15). 

C.     ELECTROMAGNETIC  SCATTERING  CROSSCORRELA  I  ION 
MEASUREMENT 

First,  a  12  inch  calibration  target  measurement  of  R  (n)2  was  performed.  To  esti- 
mate a  suitable  time  window,  a  pulse  signal  was  first  used  as  the  excitation  in  the 
anechoic  chamber.  The  measured  backscattered  signal  contains  several  constituents:  di- 
rect antenna  coupling;  chamber  clutter  and  target  scattering  which  is  modified  by  the 
antenna  responses  and  multipath.  The  target  pedestal  is  located  approximately  8  feet 
from  the  horn  antennas. 

An  additional  7  feet  of  cable  delay  line  was  linked  between  the  transmitter  power 
splitter  and  the  channel  0  sampler  to  show  the  two  signals  in  the  same  time  window.  At 
first,  a  directional  coupler  was  selected  instead  of  the  power  splitter  since  the  directional 
coupler  transmits  most  of  the  input  signal.  However,  it  was  found  that  the  directional 
coupler  provided  a  significant  attenuation  in  its  high  frequency  band.  The  result  of  the 
alignment  test  is  illustrated  in  Figure  35.  The  leading  edges  of  the  two  signals  are  aligned 
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when  the  potentiometer  reading  indicates  195.  The  time  width  of  the  sphere  response 
was  about  2  nanoseconds. 

Next,  the  noise  source  was  applied  to  measure  the  electromagnatic  scattering  cross- 
correlation  measurement.  Since  the  pulse  response  is  limited  to  about  7  GHz  [Fig.  23.], 
a  set  of  64  samples  with  a  sampling  interval  of  40  picoseconds  was  chosen  to  cover  a  2 
nanosecond  time  interval.   This  produces  a  Nyquist  frequency  of  12.7  GHz. 

Unfortunately,  the  measured  value  of  the  crosscorrelation  was  significantly  cor- 
rupted by  the  noise  and  system  errors  so  that  it  was  impossible  to  use.  Since  then,  the 
effort  was  focused  on  finding  the  source  of  the  errors  and  on  the  ways  to  avoid  it. 
However,  the  answer  is  that,  at  this  time,  we  could  not  reduce  the  error  sulficiently  to 
makes  the  estimation  useful.  In  fact,  the  light  against  the  error  noise  for  the  signal  noise 
source  spans  the  whole  history  of  this  research.  The  situation  in  the  scattering  meas- 
urement offers  signifcantly  less  output  signal  strength  than  was  observed  in  the  initial 
tests  using  lumped  filters.  This  degradation  of  signal  appears  to  have  exceeded  the  limits 
of  fidelity  needed  to  demonstrate  any  viability  of  the  measurement.  In  short,  the  results 
appeared  as  "garbage"  that  could  not  be  reproduced  on  any  two  measurements.  The 
signal  was  essentially  buried  in  the  noise. 

In  addition  to  the  description  about  the  noise  in  Chapter  II,  the  sampled  sequence 
of  the  two  channels  are  modeled  as 

Xp{n)  =  X{n)  *  pin)  +;,(//)  +  Qx(n)  +  dx{n)  +  N'x(n)  +  Dx{n)  (5. 1 7) 

Yp(n)  =  Y{n)  *  q{n)  +jy(n)  +  Qy{n)  +  dy{n)  +  sy(n)  +  uy{n)  +  Ny(n)  +  Dy(n)  [5. 1 8) 

where, 

1.  X{n)\  Input  signal 

2.  p(n)\  Impulse  response  of  the  input  channel  sampler 

3.  .;',(//);  Jitter  noise  of  the  input  channel 

4.  QJji):  Quantization  noise  of  the  input  channel 

5.  </,(//):  Drift  (noiistationary)  noise  of  the  input  channel 

6.  Nx{ri):  Other  noise  (Thermal  noise  etc.) 

7.  /),(//):  DC  bias  due  to  the  DPO  vertical  alignment 

8.  }'(//):  Output  signal 

9.  (/(//):  Impulse  response  of  the  output  channel  sampler 

10.  jy(n):  Jitter  noise  of  the  output  channel 
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11.  Qy(n):  Quantization  noise  of  the  output  channel 

12.  dy{n):  Drift  (nontationary)  noise  of  the  output  channel 

13.  sy(n):  Sampling  interval  scaling  noise  of  the  output  channel 

14.  Uy{n)\  Mispositioning  noise 

15.  Ny{n)\  Other  noise  (Thermal  noise  etc.) 

16.  Dy(n):  DC  bias  due  to  the  DPO  vertical  alignment 

It  may  be  surmised  that  the  poor  performance  is  probably  caused  from  the  nonstation- 
ary  characteristics  of  the  DPO.  as  summarized  in  Section  III.C.  This  includes  the  jitter 
noise  and  quantization  noise.  The  estimation  error  for  a  stationary  process  should  be 
reduced  by  increasing  the  number  of  sample  points.  As  shown  in  Figure  }(\  the  error- 
noise  is  not  linearly  reduced  by  increasing  the  number  of  the  samples.  Tins  reveals  the 
existence  of  other  sources  of  the  error  such  as  nonstationary  noise  (caused  by  the  drift 
problem)  which  may  give  a  sudden  level  of  noise  that  can  not  be  reduced  by  the  opera- 
tor. 

Figure  36  reveals  that  the  electromagnetic  scattering  crossco'rrelation  measurement 
for  the  12  inch  sphere  is  significantly  corrupted  by  the  error  noise  although  the  result 
came  from  the  computation  of  8192  samples.  Therefore,  it  might  be  almost  impossible 
to  perform  deconvolution  under  such  a  low  SNR  environment  unless  other  ways  can  be 
found  to  increase  the  SNR  of  the  estimated  value.  One  such  way.  which  may  solve  this 
problem,  is  the  use  of  a  new  generation  of  sampling  devices,  such  as  the  Hewlett  Packerd 
Model  54120  digitizing  oscilloscope.  This  DPO  has  a  sampling  rate  up  to  20  GHz  and 
offers  computer  controlled  time-offsets  o[  its  4  channels.  A  follow-on  investigation  will 
consider  the  application  of  the  IIP-DPO  to  this  research. 


62 


i 

i 

Transmitted  signal 
Backscattered  signal 

1 

Figure  35.      Time  Origin  Alignement  Test  of  the  8  Inch  Sphere 
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Figure  36.       The  Fluctuation  of  the  Estimated  Crosscorrelation  Due  to  the  Error 
Noise 
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VI.     CONCLUSIONS 

A.      SUMMARY 

The  main  effort  described  in  this  thesis  is  to  demonstrate  the  viability  of  impulse 
response  measurements  using  a  broadband  random  noise  source.  An  application  of  this 
system  is  radar  target  detection  and  tracking,  where  it  would  offer  the  following  potential 
advantages: 

1.  Low  cost  noise  sources  are  available. 

2.  Deception  in  radar  applications,  (jammer  or  natural  noise  appearance). 

The  first  step  of  the  research  was  to  develop  the  theory  of  noise  source  impulse  re- 
sponse measurements.  An  analog  version  of  the  crosscorrelation  was  derived  followed 
by  development  of  the  three  crosscorrelation  measurement  methods.  The  third  meas- 
urement method  (pre-subtracted  single  channel  measurement  with  arbitrary  sampling 
rate)  was  selected  after  much  effort  in  trying  to  use  method  2.  The  bias  and  variance  test 
of  the  estimated  crosscorrelation  followed.  It  was  shown  that  the  estimation  is  unbiased, 
and  consistent.  Finally,  the  effect  of  the  sampler  on  the  signal  as  well  as  the  application 
of  an  optimal  deconvolution  technique  were  considered. 

The  next  step  was  devoted  to  the  experimental  configuration.  The  characteristics 
of  the  Tektronix  DPO  and  noise  generator  were  investigated.  It  seems  that  the 
Tektronix  DPO  is  not  the  ideal  device  for  this  research  since  it  exhibits  many  problems 
for  these  measurements.  Some  of  these  are  the  drift  problem,  the  relatively  large  amount 
of  jitter  noise  (about  1/5  of  the  signal  strength),  a  sampling  interval  scaling  problem,  etc. 
An  effort  to  minimize  this  error  noise  was  attempted.  Detailed  tests  were  performed,  and 
it  was  revealed  that  the  system  crosscorrelation  included  errors  in  its  peak-to-peak  value 
of  about  1  20  of  its  maximum  signal  power. 

The  calibration  and  validation  test  were  then  performed.  An  Avantek  solid-state 
microwave  amplifier  was  selected  as  the  test  item.  The  crosscorrelation  measurements 
of  the  sampling  system  and  the  Avantek  amplifier  were  performed.  A  smoothed  fre- 
quency response  of  the  Avantek  amplifier  was  computed  using  the  Riad's  method. 
Then,  the  impulse  response  of  the  Avantek  amplifier  was  obtained  using  the  inverse 
Fourier  transform.  This  was  compared  to  results  obtained  via  a  direct  pulse  excitation 
measurement. 
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Finally,  an  electromagnetic  scattering  measurement  was  attemped  after  considering 
the  requisite  mathematical  foundation.  Four  measurements  were  needed  in  practice  al- 
though the  theory  requires  three  measurements.  The  scattering  measurement  provided 
insufficient  fidelity  due  to  hardware  limitations  of  the  Tektronix  DPO.  Major  sources 
of  error  in  the  current  system  include  nonstationary  drift,  jitter.and  mispositioning. 

B.      FUTURE  CONSIDERATIONS 

The  main  difficulties  encountered  in  this  research  resulted  from  the  use  of  a  less  than 
ideal  sampling  device  that  could  not  sample  sufficiently  fast  to  satisfy  the  Nyquist  crite- 
rion. Further,  it  has  been  revealed  that  this  sampling  device  has  many  problems  as  de- 
scribed in  the  previous  section. 

The  significant  payoff  of  this  research  mandates  that  it  be  continued  with  vigorous 
attempts  to  alleviate  the  hardware  deficiencies  encountered  in  this  initial  elfort.  The  use 
of  a  new  generation  digital  processing  oscilloscope  by  Hewlett  Packard  is  a  step  in  the 
right  direction  since  this  device  allows  computer  control  of  the  temporal  olfset  of  each 
channel  with  an  accurate  resolution  of  10  picoseconds.  The  essentially  drift- free  per- 
formance, coupled  with  automated  time-offsets  in  the  crosscorrelation  data  acquisition. 
will  allow  large  ensemble  averages  to  be  employed  with  much  shorter  measurement  times 
than  have  been  the  case  here. 
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